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(54) Property determination 

(57) A method of determining or predicting a value 
Px of a property (e.g. octane number) of a material X or 
a property of a product of a process from said material 
or yield of said process, which method comprises meas- 
uring the absorption Dj^ of said material at more than one 
wavelength in the region 600-2600 nm. comparing the 
said absorptions or a derivative thereof witti absorptions 
Djn, or derivatives thereof at the same wavelength for a 
number of standards 8 in a bank for which the said prop- 
erty or yield P is known, and choosing from the bank at 
least one standard S^, with property said standard 
having the smallest average value of the absolute differ- 
ence at each wavelength i between the absorption DjX 
(or derivative thereof) for the material and the absorption 
Djm (or derivative thereof) for the standard S^, to obtain 
P^, with averaging of said properties or yields P^ when 
more than one standard is chosen. 
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Description 

This invention relates to a method of determining or predicting by near infra red (NIR) spectroscopy properties of 
feeds or products and/or yields in physical or chemical processes or separations, in particular- involving hydrocarbons, 

5 especially in hydrocarbon refineries 

NIR spectroscopy has many advantages over other methods of analysis in refineries and can cover a large number 
of repetitive applications accurately, quickly and on line. The NIR region between 800 and 2500nm contains the totality 
of molecular information in the form of combinations and overtones from polyatomic vibrations, but f^athematical tech- 
niques are needed to exploit this information and to calculate the desired parameters. EP-A-285251 , 304232, 305090, 

'0 the disclosure of which is hereby incorporated by reference, describe the use of NIR for determining octane number, 
yields and/or properties of a product of a chemical process or separation process from analysis on the feeds to that 
process, and yields and/or properties of a product of a blending operation again from analysis on the feed thereto. 

At present, numerical methods described for modelling physicochemical properties based on NIR spectra all are of 
a correlative nature and involve relations of a regressional character between the property(ies) studied. Among these 

15 multivariable analyses are multilinear regression (MLR), Principle Component Regression (PLR), Canonic regression, 
and regression by Partial Least Squares (PLS). In all cases there is sought between the property and the NIR spectrum 
a relation which may be linear but is usually quadratic or of higher algebraic form involving regression coefficients applied 
to each absorption. The establishment of any regression requires a progressive calibration, as the approach is empirical 
and not supported by a theory, 

20 These techniques have disadvantages, the chief of which is the need for establishing a strong correlation between 
the spectrum and the property, and their difficulty in dealing with positive or negative synergy between components 
contributing to that property. For examplefor determining chemical composition e.g. LINA (linear, isoparaffln, Naphthenic, 
Aromatics) in a hydrocarbon feed to a catalyst reformer, a PLS technique based on the NIR spectra has been described 
for use. The model works well on the calibration set but the response of the models when pure hydrocarbons are added 

25 e.g. cyclohexane is not satisfactory, as the model predicts changes in isoparaffins and naphthenes the reverse of that 
found experimentally Furthermore there are other practical difficulties, mainly in the need to identify samples of families 
having the same kind of relation between the spectra and the properties to be modelled. Thus the model may be limited 
especially with a non linear relation between spectrum and property. Especially when at the edges of the available data 
the accuracy of the model diminishes. The stability of the model is also a problem, as is the need when adding new 

30 standards to do laborious revisions to give the new model, especially when adjusting to a new feedstock for a process, 
thus testing 6 properties on 4 products leaving a distillation unit requires 24 models, each of which has to be changed 
for each change of the feed not included in the calibration. 

We have discovered a new approach avoiding the above problems with correlations, and regression calculations, 
and being capable of being expanded automatically with use of a new product of different quality. 

35 The present invention provides a method of determining or predicting a value Fx, of a property of a material X or a 
property of a product of a process from said material or yield of said process, which method comprises measuring the 
absorption DjX of said material at more than one wavelength in the region 600-2600nm, comparing the said absorptions 
or a derivative thereof with absorptions Dfrn or derivatives thereof at the same wavelengths for a number of standards 
S in a bank for which the said property or yield P is known, and choosing from the bank at least one and preferably at 

40 least 2 standard with property P^. said standard having the smallest average values of the absolute values of 
the difference at each wavelength i between the absorption D|X (or derivative thereof) for the material and the absorption 
D|m (or derivative thereof) for the standard to obtain value Px, and with averaging of said properties or yields Pm, 
when more than 1 standard Sni is chosen. 

The above method can be performed without regression or correlation techniques. 

45 Thus for the performance of the method of the invention, a bank is prepared in which the NIR spectra are recorded 
at many wavelengths for a large number of standard materials, together with their properties (or those of products 
obtained by processes therefrom) determined by alternative techniques e.g. gaschromatography for chemical compo- 
sitions and yields determined by known methods. The standards are chosen to cover the area in which the method is 
to be used, so tor octane number determination, a range of gasolines is chosen of widely varying octane numbers, with 

50 different contents of lead, or other additives such as alkyi ethers and aromatics. The number of wavelengths chosen 
may be 2-1000 e.g. 5-200 or 10-20 such as 40-80 while the number of standards can be at least 100 or IOOO, or 100,000 
up to 5 million depending on property(ies) chosen. 

The wavelengths chosen may be at regular intervals such as each 1 -50 or 1 5-35nm (or each 1 -5nm or each nan- 
ometre) or may be at irregular intervals e.g. with intervals of 1-200nme.g. 1-100 or 1-50 such as 4-50 or 10-60nm, which 

55 may be gandom or chosen because of a change in the shape of the spectral curve at that wavelength e.g. a peak, trough 
or shoulder. The wavelengths_may be in the region 600-2600nm, such as 800-2600nm, in particular 1 500-2600 or 2000- 
2550nm, or 800-2000 especially 1000-1800nm for diene containing gasolines such as ones produced by cracking e.g. 
steam cracking. The wavenumbers may be in^he region 1 6,600-3840cm"i , e.g. 1 2.500 to 3840cm~i in particular 6660- 
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3840 or 5000-3900cm"i. or 12500-5000 especially 10000-5500cm-i; corresponding frequencies in Hertz can be 
obtained by multiplying this wavelength by 3x10iocm/sec. 

The absorptions for the unknown sample are compared with the absorptions at the same wavelength of the stand- 
ards, and those standards chosen having the smallest differences. The properties of those chosen standards are then 

5 averaged to determine the property of the unknown sample. The absorptions at more than one wavelength may be 
chosen, e.g. 2-1000 such as 5-100 or 10-20. 

In the method of the invention the standards chosen are those with the smallest average values of the absolute 
difference at each wavelength 1 between the absorption/optical density (or a derivative thereof) Dj^ for ttie unknown 
material and the corresponding ^sorption/optical density (or derivative thereof) Dj^ for the standard. The averages may 

10 be in respect of the mean value of Di„-D,„ (whatever its sigh i.e. absolute diference). or (Dix-Din,)2 and may be the simple 
mean vague or the differences may be weighted to take account of the different sensitivity of the absorption to the 
property at that wavelength or the different sensitivity of the spectrometer at that wavelength. For each standard in the 
bank of standards for the type of material in question, the average difference is found as described and the standard or 
standards with the smallest average differences chosen, e.g. at least 1 but preferably at least 2 such as upto 1000 

15 smallest such as 1 (or 2)-1 00 or 1 (or 2)-20 but is particular 1 (or 2)-10 and especially 2-6 smallest. Advantageously the 
average differences chosen and hence the standard (or standards) S^, chosen for the property or yield wanted are such 
that in relation to the unknown material X and each chosen standard the following functions is met 

20 < experimental error 

wherein i^^ is the proximity index and is defined by i ^ (xm) = L{D ,^-0 ^ and the experimental error is in determining 
said property or yield in the standard. The value Px of the property or yield is the same as property or yield or the 

25 average Pm if more than one standard S^, is chosen. 

In order to aid the choice of the appropriate standards, especially in relation to a large number of wavelengths for 
a complex unknown mixture, it is preferred to limit the choice to those defined by means of a minimal index. For the 
chosen standard the minimal index is at least the same as the differences between the absorptions of the unknown 
and the standards. Mathematically, this may be expressed as fiab^fiM where iM is the minimal index for the property, 

30 and iab is a measure of the deviatron (called the proximity index) at all the chosen wavelengths between absorption of 
the unknown and a chosen standard b. That measure is defined by 

i(ab)' = Ij(Di3-Dj' (I) 

35 where Dja is the optical density (or absorbence) of unknown a at wavelength i (or a derivative thereof e.g. a first, second 
or third derivative of that density), and Db is the optical density (or absorbence) of standard b at that wavelength i (or a 
derivative thereof e.g. a first, second or third derivative of that density). The value of Di is the optical density or the 
optical density difference with respect to the baseline of the spectrum at that wavelength, or the baseline interpolated 
between 2 wavelengths on either side thereof. 

40 If desired instead of the optical density Dj a normalized density W, may be used where W j =D /LD , . This normal- 
ization avoids errors due to small electronic fluctuations in the apparatus and compensates for small differences in the 
optical path between the optical cells.. In this case the proximity index is defined by 

l(ab)^ = Lj (Wi3-Wb)^ (2) 

45 

The indices can be weighted as desired for increasing resolution. One approach is to define the indices as follows. 

I(ab)'^ = lAbsvalue(Xia-Xji3)'"/ai " (3) 

so where Xj is Dj or Wj or a mathematical combination thereof, a, is the standard deviation of X for the set of samples 
considered (at that wavelength) and each of m and n which are the same or different is weighting factor which is positive 
but can be a whole number or a fraction. Other variants can be used with other weighting factors such as those involving 
the spectral experimental error e,, where ei is the reproducibility of the spectral measurement at wavelength i. The choice 
between the different options for the weighted indices may be dictated by numerical efficiency. 

55 The reproducibility of the experimental measurements in the standards may be at least 90% or 94% or 95%. The 
minimal index may be obtained from a reference standard samples set according to the following procedure, hereafter 
called the Minimal Index Procedure. The NIR spectra for 2 standard samples A and B and their property P e.g. Octane 
Number are determined. By means of eqyation (1), (2) or (3), the value of the proximity index igb is determined via the 
absorptions at a series of wavelengths, this index is applicable to the difference in properties Pg-Pb called EPab- This 
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process is repeated with other pairs o< standards c and d, e and f etc to obtain a series of Proximity Indices \^ etc vwth 
corresponding property differences EPod etc. For different values of a parameter L which is greater than the indices igb 
etc, the corresponding values of EP^b etc are averaged to give an average EPy for that value of L; the different values 
of EP i +ta/ K are then plotted on a graph against L a is the accuracy of the property determination and K is the number 
of pair's of samples for which iab is inferior to a given L. t is the Student factor at a given level of confidence. The intercept 
is then measured between the curve obtained and a line usually horizontal which is the reproducibility of the property 
level at an appropriate confidence interval e.g. 90% or more usually 95%; the abcissa portion of the intercept gives the 
minimal index i^,n, which is the minimum value of igb for which Pa=Pb within the frame of experimental error. 

From this minimal index by Procedure 1 , the standards can be chosen which have values of iZgb ^ i^min where in 
this case a is the unknown and b is a standard, as in this case the difference between Property a and Property b is less 
than or equal to cr\/2, where a is the experimental error in measuring the property Then from the property P value or 
values of the chosen standard, the property of the unknown is obtained directly or by averaging those values, usually 
the arithmetic mean, but optionally with weighting. 

The method of the invention may be used to determine more than one Property P at once, e.g. at least 2, such as 
1-30 e.g. 2-10 properties at once. Each property of the standards has a particular unweighted, minimal index, which 
may lie in the region 0-10"iOe.g. 10'2to lO'S, in particular 5x10"7 to 5x10"". if the Minimal Index chosen is the smallest 
for all the properties desired, then the same one may be used for all the properties and the standards chosen will be 
suitable for all the properties. The Minimal Index for each property may be used separately, with different numbers of 
standards chosen for each property (assuming different Minimal Indices). If desired the same Minimal Index may be 
used, which is not the smallest, resulting in some of the chosen standards (with a higher Minimal Index) giving some 
properties of high accuracy and some (with a lower Minimal Index) giving some properties of less high accuracy 

The property to be determined may be of the sample being analyzed or a product obtained from that sample e.g. 
a product of blending or cracking the sample, as the property value obtained is derived from the standards, and they 
will have been determined as needed for the eventual use. Our EP304232 and 305090 referred to above describes such 
techniques when applied to use of NIR with correlation to trending or cracking operation; the same principles apply in 
the present method. 

If the density of the standards in the data bank is sufficient to have i2ab=si2 min as is usually the case, the above 
procedure is very satisfactory But there are occasions when the bank is incomplete, because of shortage of data of 
properties in a particular area i.e. a low density of standards or the sensitivity of the property to changes in absorption 
is so small, that a very small Minimal Index is required and there may be few standards with proximity indices meeting 
it It is possible simply choose a larger Minimal Index with e.g. 1-5 times such as 1.5-2 times the Minimal Index; the 
results may be less accurate than those from a smaller minimal index. 

However, a more accurate approach with a low density of standards involves a special densification process of 
Procedure 2, in which random or semi random densif ication of the neighbourhood of the unknown is achieved by gen- 
eration of synthetic standards, based on standards already in the bank. Each new synthetic standard may be obtained 
from combinations of standards taken at random from the bank but preferably it is obtained from the other standards by 
the constraint of choosing only a mixture of N standards for which 

(Min)Cj - Uj =sCj =3(Max)Cj+Uj (4) 

and iCj = 1 (5) 

where Cj is the fraction of component j in the sample j. 

Min Cj is the minimum amount of j in the initial industrial calibration mixture or in the samples for which the method 
is to be used. 

Max Cj is the maximum amount of j in the initial industrial calibration mixture or in the samples for which the method 
is to be used. 

uj is usually between 1 and 0.01 preferably between 0.5 and 0.1 and can be fixed for each property 

The constraints over the choice of such mixtures of N standards can also be equally fixed in the spectral area from 
which the samples will be drawn in order to remain in the areas of similar chemical nature. 

The number of samples effectively drawn into the bank in this densification can be of several thousand generally 
1000-2000. The calculation time is extended without significant deterioration in the results. If no further neighbours are 
found, the trawl of new samples drawn in is enlarged. 

The spectrum of each mixture is calculated by the combination of the spectra of the standards used according to 
the formula 

S^i = LCijXS, (6) 
where Sj is the spectrum in the mixture of connponent j in the calibration matrix. 
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The properties of each mixture PMi can be calculated by a generally linear combination of the properties of the 
standards according to the formula 

PMi = £Ci,XPj (7) 

5 

where Pj is the property of component j 

In the case of non linear additive properties, appropriate mixing factors can be applied e.g. by blending factors 
or similar for density and viscosity 

Having obtained the spectrum and the properties of the synthetic mixtures, these can be used as "standards" to 
JO help determine the properties of an unknown sample in the same way as a conventional standard. 

Instead of using either of the two above approaches. 1-7. a third type Procedure 3 may be used as follows. The Q 
nearest samples to unknown X can be found from a selection from the bank samples for which the proximity index to 
the unknown sample is (V) X imm) where v is 0.1 < v < 10. (8) preferably 0.5<v<2 or Is v s5. Then by the method of 
least squares is found a generally linear combination of the standard products, which are the Q nearest samples to 
IS reproduce the spectrum of X according to the equation. 

S, = 2:CrXS, (9) 

where Cr is the coeHicient for sample R in the total Q and Sr is the spectrum of sample R. The coefficient Cr which can 
20 be normalized to Cr = 1 or not and/or optimized by the least squares route, allows an estimation of the property P^ 
according to the equation. 

P, = 5:Cf,XPr (10) 

25 where Pr is the property of sample R. 

The eventual size of the estimation error can be derived by application of Gaussian theory, also called the propagation 
error (see Eq.10). 

The above third approach can only be applied if the product X is situated Inside the maximum extensiori of the 
standard products defined by equation (8). If this is not the case. X is outside the field of the actual bank of products 

30 and escapes from the area of knowledge of the method into the area of learning. 

The densiflcatlon process described in relation to equations 4-7, or 9 or 10 Is usually applied to the method of the 
invention involving no correlation or regression techniques. However, if desired the densiflcatlon process may be applied 
to increase the number of "standards" for consideration in an NIR analytical technique involving the correlation on regres- 
sion techniques as described above e.g. MLR. The present invention also provides a method for adding an extra synthetic 

35 standard to a bank of known standards, each of which relates at least one absorption in the 600-2600nm region (or 
derivative thereof) of a known material to a known property related to that material, which method comprises choosing 
from the bank at least 2 standards for which equations 4 and 5 above are met. considering mixing the chosen standards 
in at least one proportion to produce at least one mixture for use as a synthetic standard, and estimating the spectrum 
and property of said mixture according to equation 6 and 7 respectively 

40 The spectrum and property of each "mixture" can then be added to the bank and used to develop models through 
the known correlation/regression approach, e.g. as described in the above mentioned patents. 

The method of the Invention may be applied from the spectrum of a material to determine at least one physical, 
chemical, physicochemical and/or rheological property of that material, which may be a product of a chemical or physical 
or separation process, or which may be a feed to such a process, or the method can be used to determine at least one 

45 Of said properties of a product of that process from the spectrum of at least one feed to that process, or to determine 
the yield of at least one product of that process. Each of the feed (or feeds) or products to the process may be a solid 
liquid or gas preferably at least one feed or product Is a liquid. 

Thus the method may be used for the physicochemical determination or prediction in relation to at least one feed 
or product used in or obtained by an industrial process of the refining of oil and/or in petrochemical operations. The 

50 process may be a hydrocarbon conversion or separation process, preferably a reforming or catalytic cracking or hydrot- 
reatment process or distillation or blending. In particular it may be used for determination of at least one property of a 
feed and/or the prediction and/or determination of at least one property and/or yield of product from a number of different 
processes such as processes for separating petroleum products such as atmospheric distillation vacuum distillation or 
separation by distillation, under pressure greater than atmospheric, as well as thermal or catalytic conversion, with or 

55 without partial or total hydiKjgenation. of a petroleum product, such as catalytic cracking e.g. fluid catalytic cracking 
(FCC), hydrocracking, reforming, isomerization, selective hydrogenation. viscoreduction or alkylation. 

Of particular value is the use of the method in blending operations involving the prediction and/or determination of 
at least one property of a blend of liquid hydrocarlxjns (optionally with other additives such as alkyi ethers or alcohols), 
this method including or not the determination for each constituent of the blend of a blend index for the property consid- 
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ered. tn this method as applied to blending, the fc^end indices can be obtained simply by calculation and without the 
need for preparation of standard physical mixtures other than those contained in the databank. The blend indices can 
;)e combined linearly or non linearly within the fields of stability to determine from the value of this combination a value 
for at least one property of the blend obtained. The blend may be made by mixing at least 2 of butane, hydrogenated 
steamcracked gasoline, isomerate, reformate, oxygenates (MTBE,TAME,MeOH,TBA...), FCC derived gasoline. This 
process may be repeated with numerical addition of other constituents separately to the liquid hydrocarbon base to 
determine a series of blending Indices and then determination from these indices of the properties of the multi constituent 
blend (see e.g. Ex. 2 hereafter). 

Examples of properti es that can be determi ned and/or predicted I nclude the following : for automobile fuels/gasolines, 
at least one of the Research Octane Number (RON), (\4otor Octane Number (MON) and/or their arithmetic mean, with 
or without lead additive and/or the methyl tert, butyl ether or methyl isoamyl ether and/or benzene content: 

For automobile fuels/gasolines, at least one of the vapour pressure, density, volatility, distillation curve, e.g. per- 
centage distilled at 70°C and/or lOCC, oxygen content or benzene or sulphur content, chemical composition and/or 
gum content e.g. expressed in mg/IOOmI, and/or susceptibility to lead (these properties are particularly determined for 
use in blending operations); 

For diesel fuels or gas oils, at least one of the cetane number (e.g. motor measured), cetane index, cloud point, 
"discharge point", fllterability. distillation curve, density e.g. at 15°C, flash point, viscosity e.g. at40°C, chemical compo- 
sition, sensitivity to additives and percentage of sulphur; 

For distillation products from crude oil e.g. under atmospheric pressure at least one of the density, percentage of 
sulphur, viscosity at 100°C, distillation curve, paraffin content, residual carbon content or Conradson carbon content, 
naphtha content, flash point for petrol, cloud point for gas oil e.g. light gas oil and/or viscosity at 100°C and/or sulphur 
content for atmospheric residues, and yield for at least one of the cuts, gasoline (bp 38-95°C), benzine (bp 95-1 4900) 
naphtha t»p 1 49-1 75°C, jet fuel bp 1 75-232°C, light gas oil bp 232-342°C, heavy gas oil bp 342-369°C, and atmo^heric 
residue greater than 369°C. 

For at least one of a feed or a product of a process of a catalytic cracking e.g. FCC process, at least one of the 
density percentage of sulphur, aniline point, gas oil index, gasoline index, viscosity at 100°C, refractive index at 20''C 
and/or 60^C, molecular weight, distillation temperature e.g. 50% distillation temperature, percentage of aromatic carbon, 
content of total nitrogen and factors characterizing the suitability of the feed for the cracking e.g. KUOP, crackability 
factor, cokability factor, and yield e.g. of gas, gasoline, gas oil or residue. Thus there may bedetermined the yields and/or 
properties of the different products obtained by distillation of the cracked products, such as RON and/or MON, clear or 
leaded for the gasoline cut and the viscosity at 100°C for the distillation residue. 

For at least one of a product or a feed of a catalytic reforming process, at least one of the density, distillation tem- 
perature and/or chemical composition (expressed as a percentage) of saturated linear hydrocarbon, isoparaff ins. naph- 
thenes, aromatics and olefins. 

For at least one of a product or a feed of a process of hydrogenating gasoline at least one of the density, distillation 
temperature, RON and/or MON, clear or leaded vapour pressure, volatility, chemical composition (expressed as a per- 
centage) of saturated linear hydrocarbons, isoparaffins, naphthenes, aromatics e.g. benzene, and mono/di substituted 
benzenes, olefins e.g. cyclic and non cyclic olefins, diolefins. the maleic anhydride index, and yield e.g. of at least one 
of the products obtained. 

The method of the invention may also be used with chemical reactions in which at least one product is a hydrocartwn, 
and none of the feeds or products contains an element other than carbon or hydrogen. The hydrocarbon which may be 
gaseous or liquid at 25°C. Such reactions may involve as feed or product at least one olefin or acetylene e.g. linear or 
branched, aliphatic or cycloaliphatic olefin with an internal or external ethylenic unsaturation, preferably of 2-20 carbons 
especially 2-8 carbons for alkenes or alkynes (such as ethylene, propylene, butene 1 or 2, isobutene, isopentene) or 
acetylene, and 5-8 carbons for cycloalkenes e.g. cyclohexene. The feed or product may also be an aromatic hydrocarbon 
e.g. benzene or naphthalene, optionally substituted by at least one (e.g. 1-3) alkyi or alkenyl group e.g. of 1-20 carbons, 
such as 1 -6 carbons, especially methyl, ethyl or isopropyl; examples are benzene, toluene xylene, cumene and styrene. 
The feed or product may also be a non aromatic hydrocarbon, e.g. linear or branched aliphatic or cycloaliphatic with e.g. 
1 -20 or 5-8 carbons respectively, preferably 1 -6 carbons and 6 or 7 cartjons respectively, examples are methane, ethane, 
propane, n-butane, isobutane, and cyclohexane. The feed or product may also be a diene. conjugated or unconjugated, 
aliphatic or cycloaliphatic with e.g. 4-20 carbons or 6-20 carbons respectively; examples are butadiene and isoprene 
and cyclohexadiene. Examples of the reactions are hydrogenation (e.g. butadiene to butene-l or 2 or cyclohexene to 
cyclohexane) dehydrogenation (e.g. ethane to ethylene or ethyl benzene to styrene), isomerisation (e.g. butene- 1 or -2 
to isobutene. or pentene-1 to isopentene) alkylation (e.g. benzene with ethylene to form ethylbenzene and/or styrene, 
or isobutene with butane to form iso octane), and cracking. 

In addition to the use in getrochemical operations, the method is of wider application and may be applied in the 
pharmaceutical industry such as the production of pharmaceutically active compounds for use as medicines e.g. by 
fermentation, and in the perfumery industry for making perfumes and fragances, especially in their blending and control 
thereof. The method may also be used in the feod industry e.g. in brewing to control fermentaion processes, in fermen- 
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tation to make wine and quality control thereof, and control of food production e.g. sugar and water content in fruit juice 
and in control of maturing processes for fruits and vegetables In each case the method may be applied to determine a 
property of the sample tested or product from that sample e.g. a fermentation or blended product preferably on line and 
especially with continuous feed back from the results to control the production process. 

5 In each of the above processes the property or yield of a product determined or predicted by the method of the 

invention can be compared to the desired figure and notice taken of any deviations by adjusting the parameters of the 
process e.g. proportion or nature of feed(s) and/or temperature/pressure etc to bring the property back to the desired 
figure. This control of the process, which may be a blending, separation or chemical process, is usually performed with 
a micro computer which is linked to the spectrometer and also performs the search for the standards Sm. The inline 

10 control of the process is very efficient and very fast. 

The present invention also provides an apparatus suitable for carrying out the method of the invention comprising 
an infra red spectrometer and a computer wherein the infra red spectrometer is linked to the computer programmed in 
such manner that the property or yield may be determined continuously and in real time. The spectrometer is suitable 
for measuring spectra in the 600-2600nm wavelength range and can be linked to a signal processing device to allow 

75 numerical treatment of the spectrum, preferably by Fourier Transformation. The spectrometer receives at least one signal 
from a vessel containing product or from a feed or product line. The information obtained can be used as an information 
vector for the computer which is programmed to determine the property or yield e.g. via calculations on the proximity 
indices in relation to standards. Conveniently in relation to a process, the computer may be used in a closed loop feedback 
control system for controlling processing equipment e.g. changing the process parameters in response to variations in 

20 the property and/or yield of product from the desired value, from measurement of more than one absorptions in the NIR 
spectrum of the product and/or feed. 

The benefits of invention allow improvements in modelling with the following areas, identification and classification 
of novel products, simultaneous estimation of all of P properties on a sample without the need for generating P different 
models, and with the option of automatic upgrading of the model, the method being self learning or adjusting. The method 

25 of the invention overcomes the difficulties with the classical regressional approach, in particular avoiding all difficulties 
with numerical stability of the models, allowing easy and rapid identification and classification of a sample of a product 
analyzed by spectral recognition and then instant conclusions as to whether the sample is known or unknown, allowing 
simultaneous determination of many properties and whether the property is simply additive or synergetic in relation to 
a blend composition: the latter is particularly useful for different blend indices and the indices considered. 

30 The method also allows an extension of the field of application of the method without the need to rewrite the model , 
apart from the need to integrate the new samples which are inside or outside the previous field of validity of the method. 
This possibility of automatic learning, which is not possessed by traditional regression techniques, is a decisive advan- 
tage in the framework of continuous inline industrial control processes, because it allows the return of the industrial plant 
operations to the model in a certain and rapid manner in a minimum time and with all the properties considered in the 

35 model. In contrast classical regression methods would necessitate the redevelopment of all the models, which is long 
and laborious without being able to guarantee the result of the new mode! obtained, because a new validation period is 
necessary; in addition during the redevelopment of the model any commercial refinery use of the model is very limited. 
Furthermore, the method of invention allows equally the easy extension to a number of properties, which are simply 
incorporated into the known bank. 

40 This remarkable possibility is true not only for conventional properties such as physical chemical and/or rheological 
properties, but also for complex ones (such as octane number). Also it is possible to quantify by the process the response 
or susceptibility to lead of automobile fuels as well as the response to additives such as nitrates, of fuels used in diesel 
engines. The methods of the invention equally allow application of the models from one apparatus to another and from 
one spectral region to another, where conventional regressive method cannot give satisfactory solutions. This apparatus 

45 portability is made possible by the fact that the differences between diHerent spectra are the same in one apparatus as 
another, for the same type of spectrometer being considered (e g. network scatter, Fourier transform, accousto optical 
system AOTS, diode array etc). This portability between spectral regions depends on the fact that as the spectral regions 
are intercorrelated, the relations between the spectra are maintained between one another 

The invention is illustrated in the following Examples in which the Minimal Index is calculated according to the Minimal 

so Index Procedure described above. Mathematically the steps concerned are as follows. 

For each couple of standard samples i, j, the Proximity Index iij is determined from the NIR spectra by use of equation 
1 , 2. or 3 and the properties are measured. For each Proximity Index is calculated the absolute difference EPy between 
the laroperties of the samples. The Minimal Index for property P is obtained from the average (EMpL) of EPy for different 
values of L when L gij. Thus the EM pL = 1/K i: z: EP y for each of K samples for which ij sL. 
55 EMp(L)+ta(M) is plotted against the proximity index and in addition there is plotted the reproducibility of the standard 
method at a given level of^ponfidence, as defined in the Minimal Index Procedure above. The intercept of the curve from 
EMpL and the reproducibility give the upper limit i.e. the Minimal Index. 
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For the Examples the data is expressed in Tables in a form as shown below in which the data is as follows. 



10 



15 







Absorption 






Weighting 


Unknown 


Estimated 


Standard A 


Standard B 


Proximity Index 












Wavelenc 


}th X 




































































Property I 














Property j 














Property m 















The wavelengths chosen are shown in columns 1 and 2. 

Column 3 gives the weight loading associated with each wavelength for the proximity index for the standards; 1 
denotes no loading. 

25 Column 4 shows for the unknown sample the absorption at the various wavelengths and at the bottom the properties 
of that sample determined by standard methods. 

Column 5 shows for the unknown sample the estimated values of the properties and the absorptions using the 
method of the invention based on the properties and absorptions of the chosen standards. 

Columns 6, 7 etc show the values of the absorptions and properties for the standards chosen from the bank. Line 
30 2 give the value of the proximity index between the unknown sample and each of the chosen standards. 

Example 1 

Determination of Octane Number and other Properties of a Motor Fuel 

35 

The NIR spectra between 4800 and 4000cm~i of a superfuel ID and a number of standard superfuels of known 
properties were measured. The base line was taken at 4780cm"i though similar results would be obtained with baseline 
drawn between 2 or more points. The absorbances were normalized . 

By the Minimal Index Procedure described above, with use of equation 2 and non weighting of the absorbences the 
40 Minimal Index (Ml) was calculated to be 1 x 10"4. Following reference to the bank of data on superfuels and use of 
Procedure 1 . 3 standard samples were found with a proximity index with respect to the superfuel of less than Ml . The 
properties of these standards are shown in Table 1 . From the properties of the standard samples, octane numbers (RON 
and MON), vapour pressure (hpa) volatility, percentage distilled at 70°C and at 100°C, gum content (in mg/ml), and 
content of sulphur, benzene (vol7o) and MTBE were calculated for the superfuel by taking the arithmetic mean of the 
45 values for the 3 chosen standards. The estimated results are compared with the measured results. 

All the properties were obtained from the single NIR measurement on the unknown superfuel and without any 
regression calculations, and with an accuracy in agreement with the reproducibilities of the reference methods. Other 
properties can be determined in a similar way. 

50 Example 2 

(a) Production of an Unleaded Mixed Fuel from 6 Components 

A target SUPER98 superfuel of the properties given in column 3 of Table 2a1, was to be obtained by mixing the 
55 remains pf a tank of finished gasoline with 5 components, butane, hydrogenated steamcracked gasoline HEN, isomerate 
ISOM, reformate (REF) and MTBE. NIR absorptions at 4800-4000cm"i measured with a Fourier Transform spectrometer 
were measured, with a base iTne taken at 4780cm"i and absorbances normalized. Results are in Table 2a1 . 

Mathematic calculations were done with a computer to mix the spectra and properties of the 6 components to 
reproduce a finished product. 
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5% MTBE (on target fuel) (i.e. 4-76% in the final mixture) was "added" mathematically to a spectrum of the target 
fuel to give a mixture whose NIR spectrum was noted. The Minimal index was 1 x 10"'* determined as described above 
from the finished gasoline. 3 standards 2A, 2B and 2C were found with proximity indices with respect to the mixture, 
without weighting, and hence by averaging the properties of the standards the properties of. the mixture were obtained. 
5 Table 2a.2 shows the spectrum of the mixture, the 3 standards and the estimation for the mixture as well as the properties 
of the standards and the estimated figures) . The process was repeated with addition of each of the other 4 components 
to the spectrum of gasoline target. 

On the basis of the figures obtained, the blending index for each property was found according to the linear formula 

,0 I " (mix) = [(1+a)xP(mix)-P(ref)]/a 

where P (mix) is the blending index for the ingredient in the mixture in relation to property P 
a is the percentage of ingredient in the mixture 

P (mix) is the property of the mixture (ingredient + gasoline) added) estimated by the process. 
75 P (ref) is the property of the reference target gasoline. 

The blending index for addition of N/ITBE is shown in Table 2a3. 

In order to obey the linearity law here, it is necessary to limit the additions to not more than a quarter of the minimum 
to maximum range of the constituent studied in the industrial mixtures. However for concentration less than 20% such 
as for these oxygenated compounds, addition of 5% is acceptable. 

20 The process with MTBE added to the gasoline was repeated with the other 4 components (and on the basis of 
linearity in the blending as with MTBE) to obtain blending indices for them as well (see Table 2a3). Then with the blending 
indices for each property for each ingredient, one can calculate the relative volume fractions needed to give the desired 
properties for the Superfuel 98 and hence the blending order. The 6 components were then mixed in the desired pro- 
portions and then properties of the mixture tested and compared to those estimated by the method of the invention from 

25 the components present (see Results in Table 2a 4). In the estimation of the products and the comparison with the bank 
of standards, the Minimum Index was 1 x 10"". 3 standards 2D, 2E, 2F were found with suitable proximity indices from 
which the properties of the superfuel were estimated by averaging as described in Procedure 1 There was good agree- 
ment between the properties obtained via the blending order, these measured on the fuel made and those estimated 
by the method of the invention. The differences are very small and in the area of reproducibility of the standard methods. 

so 

Example 2b 

Production of a 5 com ponent leaded Superfuel mixture 

35 A target superfuel of the SUPER 97 type had with 0. 1 5g/l of lead tetraethyl and having as specification an RON of 
97, an NIR spectrum as in Col 3 of Table 2bl below and other properties as given in col 3 of Table 2b1 below. There 
were available 4 components (HEN, 150M, REF and an FCC cat cracker gasoline) and the remains of a tank of finished 
refined gasoline for making the target fuel. The NIR spectra of these 5 components were measured as in Ex 2. The 
results are in Table 2b1 . 

40 As in Ex.2a, mathematical calculations were done with a computer to obtain the spectra and properties of 5 com- 
ponents to reproduce a finished product. Proximity indices with respect to standard samples were calculated based on 
normalized absorbencies which were not weighted. The method of the invention was used to find appropriate standards 
, using the procedure of artificial mixtures as described in Procedure 3 and equation 8 above in which v was 1 and with 
a Min. Index of 2 x 10"4, the latter having been calculated for standard fuel mixtures as described above. Table 2b2 

45 describes the results of addition of 5% of the FCC gasoline to a reference Super 97 gasoline targetas well as the 3 
standards 2G, 2H, 2J found by the method of this invention, from which the estimated properties were found. The same 
procedure was performed with the other components. 

The blending indices were found in the same way as for Ex 2a. with the results for FCC gasoline in Table 2b3 and 
for the other components in the same way. The spectral blending index (for the linear area) is obtained for each property 

50 as shown in Table 2b3. A blending order was also calculated, as in Ex 2a, the results being in Table 2b3^ 

The process of Ex 2a was repeated but with the alx)ve components and a Minimal Index of 2 x 10 ^. The results 
are in Table 2b4. 

3 standards 2K, 2L, 2M were found with appropriate proximity indices, which allowed the properties of the product to be 
estimated by averaging. Again good agreements is seen between the properties estimated from the blending order and 
55 those measured on the product made, and also between the same properties measured and those measured by the 
process. The differences seen are very small and in the area of reproducibility of the standard methods. Other properties 
can be obtained in a similar way. 
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Determination of cetane index and other properties of a gas oil 

5 The properties of an unknown gas oil 3A were desired. The method of this invention was applied with respect to a 

bank of known standard gas oils with known NIR spectra. The NIR spectra were obtained by F T spectrometer in the 
4800 - 4000 cm"i region {with 4780 cm"i baseline and were normalized] The proximity indices were calculated on the 
basis of Equation 2, and the Minimal Index was 2.5x1 0'S (estimated from standard gas oil data as described above). 
The bank of standards was sufficiently dense for there to be found 2 standards 3B and 3C inside the sphere with proximity 

10 index less than 2.5 x 10"6. Table 3.1 gives the details of the spectra and properties of the unknown oil A, and the 
standards and the estimated spectrum and properties, obtained by averaging. All the properties were obtained with an 
accuracy in agreement with the limits of reproducibility of the reference methods. Other properties can be obtained in 
a similar way. 

f5 Example 4 

On line prediction, based on NIR spectra on a mixture of crude oils fed to an atmosphere distillation unit, of yields 
and properties of the different distillation cuts such as gasoline {38-95'C) benzine (95-149°C) naphtha (149-175°C, jet 
fuel (175-232''C) light gas oil (232-242°C) heavy gas oil (342-369°C) and atmospheric residue (bp) .369''C). 

20 An atmospheric distillation unit in a refinery was fed with a charge 4C which was a mixture in w1% of the following 
crudes, RUMASHKINO 81%, Iranian Heavy 18%, Iranian light i%. 

Yields of various "distillation cuts were desired, the boiling ranges being given above, as well as key properties of 
each cut as described in Table 4. 1 , NIR spectra were measured as in Ex 1 on the crude oil. Min. Index was determined 
from NIR spectra on standard crude oil (as described above) and was 2.6x1 0"6. The method of the invention was applied 

25 using Procedure 3 and equation 8, in which v was 1 . to the bank which was sufficiently dense for 2 standards 4A and 
4B to be found with small enough proximity indices. These standards contained (wt%) (for 4A) Romashkino 52% Iranian 
Heavy 29%. Arabian Heavy 1 1%, Kuwait 4%. Arabian light 2% and Iranian light 2%) and (for 4B) Iranian Heavy 78%. 
Romashkino 21% and Arabian Heavy 1%. The data in Table 4.1 shows the observed properties as well as the yields of 
the cuts and their properties. The results obtained by this procedure were extremely satisfactory, the differences observed 

30 being in accordance with standard methods of measurement. Other properties can be obtained in a similar way. 

The yields and properties of the distillation cuts remarkably were obtained directly on the basis of the NIR spectra 
of the feed and in line without regressional type calculations. 

Example 5 

35 

Determination in line of the properties of a mixture of crude oils 

Other properties of the charge mixture of crude oils of Ex 4 were sought, based on the NIR spectra determined as 
in Ex 4. The method of the invention was applied as in Ex 1 with the Minimal Index in all cases being 2.8 x 10"8. Two 
40 standard aude petrols 4 A and 4B were found in the bank by using Equation 2. The results are shown in Table 5.1 . Other 
properties can be obtained in a similar way. 

Here too the method demonstrates its capacity to predict all types of properties without any regression type of 
calculation requiring fastidious calculations. The results generally, as in the other Example, were in accordance with the 
results obtained by the reference methods, the deviations being found in the limits of reproducibility of the same methods. 

45 

Example 6 

Determination of the Properties of a feed to a reformer 

so A feed 60 to a reformer unit was analysed by the method of the invention as described in Example 1 with the NIR 
spectra recorded at 2000-2500 tim. the absorbancies normalised and not weighted. The NIR spectrum was compared 
by the method of Procedure 3 and equation 8 (wherein v is 1 ) with a Minimum Indexof 2 x 1 0"4, which had been previously 
calculated as described abvoe from NIR spectra on standard reformer feeds. Three standards 6A. 68 and 6C from the 
reference feed bank were found with small enough proximity indices; details of the spectra of the feed and the standards 

55 are giveq in Table 6. 1 , together with 5 properties estimated for the feed by averaging the corresponding values of those 
standards. The actual properties of the feed were measured for comparison; the measurements were by traditional 
methods (gas chromatography and density), the former necessitating laboratory determination for several hours, com- 
pared to the present NIR process which gave the same results in a few minutes and on line (real time in the unit) and 
with better reproducibility. 
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The process allows the obtaining of a result with remarkable economy while avoiding having to produce 5 regressive 
models. The differences between the 5 properties as estimated and as measured experimentally are in agreement with 
the reproducibility of the known reference methods, namely 1.5% for gas chromatography for chemical compositions 
and 2%o for density The method can be equally applied for other properties such as ASTM distillation temperature 
5 curve for the feed. 

Example 7 

Determination of the properties of a feed to an FCC unit, as well as the yield and properties of the products obtained 

w 

The NIR spectrum of the alxjve feed 7D was measured at 4800-4000 m'l , with base line at 4780 cm"i , normalisation 
of the spectrum and no weighting. The procedure 3 was used with equation 8, with v = 1 , and the M\n. Index of 2.5 x 
10'6 the latter having been previously calculated as described above from NIR spectra on standard FCC feeds of known 
properties. 

J5 The properties of the feed charge 7D sought were listed in Table 7. 1 and included factors characterising the charge 
to the FCC unit, such as KUOP, crackability and cokability. The KUOP or Watson factor is defined as 

KUOP = ^Ve/density 60/60 

20 where 0 is boiling point on a Rankin scale (Absolute Fahrenheit scale) and density 60/60 is the density of the feed at 
eCF compared to that of water at 60°F 

The cracking unit operated under the following conditions; riser inlet temperature 250°C, riser outlet temperature 
525°C, MHSV (Mass Hourly Space Velocity) 78 kg/h per kg, C/O ratio 6.6, activity of catalyst 65 (in Microactivity Test). 

The cracking gave a gasoline cut defined by ASTM distillation with initial point of 38'C and 90% distilled at ISO'C 
25 and a residue defined by ASTM distillation with 10% distilling at 385°G. 

By application of Procedure 3 to the bank of samples of FCC feeds 2 standards were found namely 7A, 7B and the 
properties and yields estimated as shown in Table 7.1. The results were all in line with the accuracy based on the 
reference methods, as well as in line with the properties and yields actually meansured. Other properties of the charge 
or products can be estimated in a similar way. 

30 



On line determination of properties of the feed to a gasoline Hvdroaenation unit. 

35 The gasolines obtained from steam cracking units have the inconvenience of containing non negligible amounts of 
unsaturated dienic compounds, which have the effect of inducing and encouraging formation of gums which are unde- 
sirable in motor fuel. These gasolines are therefore selectively hydrogenated to eliminate the dienes without at the same 
time hydrogenating other unsaturated compounds present in the gasoline such as monoolefins and aronriatics. The 
control over these dienes is therefore essential not only for the final quality of the fuel {principly RON and MON) but also 

40 for the hydrogen consumption of the hydrogenation unit. 

Units for Hydrogenating gasolines from steamcrackers are generally coupled to a downstream distillation unit to 
separate a gasoline from a light cut (95% distillation by about 75°C) and one from a heavy cut (initial point about 95°C). 
before-extraction of the benzene in the core cut and recycle of the extraction residue from that cut called raffinate. 
It was desired to determine by the process of the invention the properties of the gasoline from the steam cracker. 

45 which was a feed to a gasoline hydrogenation unit. NIR spectra were obtained on the feed on line at 1000-1600 nm 
using a scatter dispersion spectrometer The absorbences were normallised, but the data was not weighted for use in 
Equation 8, in which v was 1 and 'min was 2.5 x ^0'5 (the latter having been determined from NIR spectra on similar 
feeds of known properties). 5 standards 8A-8E were found in the search using Procedure 3, and the properties of the 
feed calculated therefrom by averaging were all in agreement with the measured properties of the feed. The results are 

so shown in Table 8.1. 

In addition the chemical composition of the feed was obtained with great particularity allowing a distinction to be 
made for example between cyclic and non cyclic olefins as well as benzene and mono and di substituted aromatics. 
Equally by the process potential yields were obtained of the distillation cuts after the selective hydrogenation of the 
gasoline. All the properties were obtained with great accuracy within the limits of the experimental reproducibility for that 
55 kind (Of property. 

Other properties can (Redetermined such as Octane Indices for the different cuts or temperatures of ASTM distillation 
curves for the gasoline. 
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Example 9 

Method tor use when the density of standards in the bank is insufficient 

5 The MON level for a reformate 9A was sought. The NIR spectrum was measured at 4800-4000 cm"i with a base 
line at 4780 cm"i : the spectra were normalised. With reference to NIR spectra on reformates of known properties the 
IWinimal Index was found by calculation as described above to be 2x1 G"5. The proximity indices of reformate 9A and 
known standards were determined by Procedure 1 . The results were as given in Table 9. 1 . 5 standards 9B-F were found 
from the reformate bank with proximity indices low in relation to the reformate 9A, but insufficiently low to be less than 

TO Minimal Index, as the density of the iDank was too small. It was thus not possible to calculate the properties with the 
accuracy desired. Procedure 1 using Equation 1 was replaced by Procedure 2 using Equations 4-7, with in Equation 4 
values of Cj between -0.3 and +1.3, in order to increase the density of "standards" in the bank by providing new synthetic 
mixtures. 

Tables 9.2 and 9.3 show the results obtained, showing in Column 3 the absorbancies and properties for the "stand- 
15 ards" (MCI , MC2) obtained by this densif ication, and with small enough proximity indices. Col.4 and subsequent columns 
give the absorbances properties of the standards 9B, 9D and 9G in the reformate bank used to generate the new "stand- 
ards". Line 2 in these Tables show for each standard the fraction retained in the mixture to generate the new "standards". 
This fraction can be negative, but comprises between -0.3 (or -0.4) and -hi. 3 (in Eq.4). 

Using the data on MCI and MC2 as "standards", the properties of the reformate 9A were calculated by averaging 
20 (as shown in Table 9.4). The calculated MON of reformate 9A accords well with the experimentally measured figure, 
and is inside the limits of reproducibility of the standard method. The process can be used in a similar way for other 
properties. 

The method of the invention equally allows immediate automatic upgrading of the bank by automatic integration of 
the new samples. The process with the proximity indices allows consideration as a standard of all the novel "standards" 

25 introduced into the bank. This property is remarkable because it allows very rapid determination of properties in the 
case of novel products not recognised in the bank and then the gaining of precious time in the adjustment of operating 
conditions for the manufacturing unit. 

Table 9.5 shows that a novel "sample" measured immediately after incorporation of the above unrecognised sample 
9A in the databank, now used as a standard, is recognised and is perfectly calculated for the totality of its properties 

30 and without any modification nor intervention on the used models. It is important to note the superiority of the procedure 
over classical regressional models. The latter are incapable of predicting properties of samples not included within their 
application range or predict them with a non acceptable error, and therefore would need to be reactivated by the necessity 
to remake the model (one for each property) and this without guarantee of success, and with the commercial plant 

35 
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functioning blind during the recalibration period. 



TABLE l.I 

Dcterrainntion of Octiine Indices and other Properties in an autbniohile fuel 




^trBE 



TABLE 2a.l 

NTR Spectra of Unleaded mixed fuel and base fuel and additives 
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TABLE 2a.2 
Effect of addition of MTBE on the Super Fuel 



TABLE 2a.3 



Blending Indices and Blending Order 




(blending order) 


Base Fuel 


Butane 


HEN 


ISOM 


MTBE 


REF 


Volume Fraction 




19.30% 


4,10% 


31.70% 


32,10% 


5,60% 


7.2% 


















RON clear 


99.4 


100,2 


97,5 


103,6 


93,0 


109,6 


100,7 


MON dear 


88,2 


88,0 


88,5 


88,0 


86.7 


98.7 


88,42 


Vapour Pressure 


709,1 


767,0 


4700,0 


98,7 


923,0 


460 


208,2 


Volatility 


972,4 


975,0 


5000,0 


212,0 


1430,0 


718.3 


177.3 


%Dist 100°c 


54,9 


50,2 


200.0 


-5,0 


97,5 


99.7 


24 


%Dist 70°c 


37.8 


31.3 


142.8 


-12,9 


84,2 


63.8 


-8,3 
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TABLE 2a.4 

Comparison of the result obtained via the blending order and those of the product 
obtained 



30 



TABLE 2b.l 

NiVIR spectrum of Target Super Fuel, and base stocks available 

35 ^_ 









Gasoline 








REF 


































41Z0 




0, 0014234 








0,0045838 


3,0004093 


0, 001657 






0, 00S2323 


0, OOS-0594 






0, 037011 


D,002C427 


0, 024445 




2155 












3,0021545 




4615 




0,013236 








0,046361 




0,032231 


4585 


2181 


0, 0058247 






0077879 


0,028857 




0,023754 








0,010375 




, 014512 


0,015436 




0,012233 




2242 




0,014075 


0, 01S828 






0, 017515 


4385 






0,093263 


0, 095498 


0,078838 


0, 095201 


0, 085595 


4332 


2308 




0, 1403 


0,1543 


0,092697 


0, 14681 




4305 








0, 13213 


0, 051636 




0, 11073 


AZiO 


2347 


0, 11073 


0, 11062 


0, 11559 


0,086428 


0, 11763 


0, 093542 


4210 


2375 


0,068913 


0,065012 


0, 073406 


0, 052861 


0, 072136 






23S8 


0,069683 




0,073798 


0. 049022 




0, 05775 






0,071227 


0,071736 


6,057143 


0, 045477 


0,085471 


n, 05255 


4105 


2436 


0,070003 


0,07081.9 


0,065C03 




0, 082203 


0, 056175 




2463 




0, 10186 


0,089357 






0, 1243 




2475 








0, 1296 






































MCN cloar 
















Vapour Pressure 
















Volatility 
































•-Dist 70*c 
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TABLE 2b.2 

Effect of addition of 5% FCC gasoline on Super 97 Product 




In this Table 3.41E-05 means 3.41x10-5 



30 



TABLE 2b.3 



35 



45 



Blending Indices and Blending order 




(blending order) 


Gasoline 


FCC 


HEN 


ISOM 


REF 


Volume Fraction 




12,93% 


31,85% 


8,46% 


37,55% 


9,20% 
















RON clear 


97.7 


97.9 


91.6 


101.1 


101.5 


100. 


MON clear 


86.3 


86.4 


79.9 


83.5 


92 


87.9 


Vapour Pressure 


589.9 


6480 


388.1 


137.7 


930.0 


235.0 


Volatility 


952.7 


9685 


628.2 


274.2 


1559.0 


204.7 


%Dist 100°c 


65.5 


63.2 


57.5 


6.3 


96.0 


26.4 


%Dist IQX 


45.4 


43.1 


30.9 


-10.7 


83.3 


-4.3 



50 
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TABLE 2b.4 

Comparison between the results from the blending order and the product obtained 



TABLE 3.1 

Determination of cetane index and other properties of a gas oil 









Gas Oil A 












Measuced 


Estimated 








PrcximiCy 






1, 71E-C6 


1, 39E-G5 


















(cm-l) 
















2118, 6 




0,000120383 


0,000126618 


0,000139825 
















0, 002011949 




2155, 2 




0,003434147 


0,003415109 


0,003438675 


0,003391543 








0,004544314 


0,004490799 


0,004476561 


0,004505037 




:iEl, 5 




0,004729896 


0, 0046754 


0,00463465 


0,004715149 








0,007119883 


0,006932337 


0,006908771 


0, 006955903 




2242,2 




0,010349409 


0,010133388 


0, 010064653 


0,010202122 










0,074925207 




0, 074920298 








0, 158677852 






0, 157499551 


4304 






0,101824835 


0, 102266697 


0,102217602 


0, 102315793 


4260 


2347, 4 




0, 131871507 


0,131331453 


0, 1313B0548 


0,131282357 


4208 


237S, 4 




0, 08862786S 


0,088637684 


0,08855913 


0, 088715237 


4168 


2399, 2 




0,092899205 




0,093105408 


0,092913844 


4132 


2420, 1 




0, 084503? 12 


0,08475911 


0,084827844 


0,084690376 


4104 


2436, 6 




0,081106377 


0,081297851 


0, 081361676 


0,081234027 




2463. 1 




0, 08642837 


0,086487285 


0, 086487285 


0,086487285 




2475,2 




0, 067192608 


0,05740372 


0,067413449 


0,067388991 


c«cane Index 




52 




53, 3 


51, 7 










52, 1 


53, 4 


Density IS'C 








0, 3385 


0,8432 








57,5 




55 










0, 23 


0,27 


Cloud Point 






5,5 


5 


6 


Filtrability 






0,5 


1 


0 






3, 1 






3, 7 
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TABLE 4.1 

Determination of yields and properties of cuts from distillatioa of mixture of crude feed 
oils 



20 



30 























































































21S5.: 












4616 












































2242,2 












438 4 














4332 














43C4 


2323,4 




Oi 104631107 








4260 


2347, 4 


















0, 067615393 


0, 037751054 


0,087838988 


















413i 


2420,1 




0, 064648925 


003470632$ 
















0,082389016 


0, 08 23 409s: 








0,087068028 




0, 087475 




4040 






0,067784043 


0,056366138 


0,068439449 
















'4 Gisoline 












% Beniine 










7,2 














« Petrol 












i light eras oil 1.30 




15, 9 








% Heavy gas oil 












» Residue BAT 












% Paraffiaes Mapbta 






52, 1 


















Cloud poi:it ISO 




-3,1 






-7, 5 


« Sulphui RAT 












Viscosity 100 "C SAT 




53,16 




48, 53 
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40 
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TABLE 5.1 

On line determination of properties of a mixture of crude oils 































Proximity 








































467: 


:i40,4 






0, 001747949 




0, 00172484 












oooa^sL?" 












0, 003826312 








2181, 5 




0, 003794844 


0, 0037S6319 














0, 006142161 






4460 


2242, 2 




C, 009258476 


0,009152272 


0,009273227 


0,009031317 


4384 






0,073089814 


0,077868554 


0,077637461 


0,078099647 


4332 




308, 4 






0, 15773336 


0, 15773336 


0, 15773336 










0, 104631107 


0, 1044836 


0, 104139419 


0, 104B27792 


4260 




347, 4 




0, 130690546 


0, 1301S8858 


0, 13039S533 


0, 130002183 


4Z0S 




376, 4 




0,087915393 


0,087717055 


0,087805559 


0,087628552 










0,091208037 


0, 090844188 


0, 0S0844198 


0, 090844188 










0, 084648925 


0,08467351 


0,084619424 


0, 084727595 


4104 








0, 0824855 




0,082357651 


0, 082308493 


4060 












0,08744171 


0, 087648219 






475,2 






0,06833965 


0,068413403 


0, 068265897 






0, 8663 


0,86555 


0,8646 


0, 8665 


Sulptur 




1,6 


1, 65 






Viscosity lOO'C 




2,27 


2,265 


2, 36 


2, 17 


4 Conradson Carton 




4, B 


5 


4,8 


5, 1 






5 




4,9 


5 
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TABLE 6.1 



Determination of the Properties of a reformer feed 




Weight 


Feed 6D 


6A 


6B 


6C 






Measured 


Estimated 








Proximity Index 








5.7E-4 


8.2E-5 


9,7E-5 
















A(nm) 














2210 




0,04624 


0,04659897 


0,04648 


0,04671 


0,04659 


2260 




0,18118 


0,18154437 


0,18233 


0,18132 


0,18085 


2266 




0,25391 


0,25482278 


0.25605 


0,25439 


0,25386 


2276 




0,33866 


0,33942652 


0.3412 


0,33844 


0,33857 


2286 




0,33776 


0,33747772 


0.3395 


0,33634 


0,33652 


2307 




0,54602 


0,54558172 


0.54286 


0,54922 


0.54375 


2328 




0,38819 


0,38770261 


0.38812 


0.38791 


0.38685 


2344 




0,4557 


0,4561672 


0,45568 


0,45592 


0,45717 


2376 




0,31751 


0,31727184 


0.31483 


0.32067 


0,31543 


2397 




0,33674 


0,33644352 


0.33466 


0.33858 


0.33561 


2408 




0,31787 


0,31746329 


0,31737 


0,31904 


0.31525 


2418 




0,32524 


0,32334235 


0,3244 


0,32341 


0,32186 


2437 




0,34758 


0,34790932 


0,34915 


0,34543 


0,34996 


2457 




0,38142 


0,38057046 


0 3793 


0,38076 


0,38195 


% Linear Saturated 


33.0 


32,6 


32,5 


32,4 


32,9 


% Isoparaffins 


30,1 


30,8 


31,4 


31,1 


29,9 


% Naphthenes 


29,3 


29,2 


29,2 


28,5 


30 


% Aromatics 


7,6 


7.4 


6,9 


8 


7,2 


Density 


0.7151 


0.7158 


0.7152 


0,7167 


0.7155 
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TABLE 7.1 

Determination of Properties of feed to FCC reactor and yields .1nd properties of 



75 





















Measured 


Estimated 








Proxuaity 






1, lOE-06 


l,28E-06 


1, 30E-06 
































;il3, S 




0, 00024017 


0, 000283004 


0,000236346 


0,000327;;; 








0, 002238801 




0,001890879 


0,00212585 




;i55,2 




0, 004237234 


0, 003903227 


0,003874117 


0,00393233b 


4612 


2168,3 










0, 005079102 


4584 


2181, 5 




0, 005332797 


0, 005055095 


0,005031089 


0,005079102 








0, 007970887 


0, 007756354 


0,007744263 


0,007768446 




2242,2 




0, 011303264 


0, 011210967 


0,011212199 


0,011209736 


4384 


2281 




0,072994455 


0,07273491 


0, 07292398 


0,072545839 




2308, 4 




0, 1520S7S43 




0,151945649 


0,152373047 


4304 


2323,4 






0, 100397569 


0,100601923 


0, 100193214 


4260 


2347,4 




0, 131209247 


0, 131514201 


0,131487607 


0,131540794 


4212 


2374,2 




0, 091618024 




0,091564633 


0, 091681751 




2399,2 








0,09427582 


0, 094370104 










0, 086675314 


0,086678538 


0, 08657205 




243a, 6 






0, 081916022 


0,081581133 


0,081850912 








0,094267922 




0,084444043 


0,08419206 


4040 


247S,2 




0,06911082 


0, 05914S752 


0,069239547 


0, 069053957 






0,926 


0, 9225 


0,922 


0,523 


i Sulphui 




1, 97 






1, 87 


Aniline Point 




83, 5 




78,2 


88,2 






E,3 




8,7 


9,5 


Temp. 50% dxsCilled 






464 




471 


KUOP 




11,8 


11,85 


11, 3 




Mol. Weight 




450, 6 


449,95 


434, 5 




i Axomatic carbon 






21,2 


21, 6 




CRACKS3ILITY 




2, 47 


2,57 










1, 01 






1, 01 


GASOIL INDEX 




1, 55 


1,515 


1,54 


1, 45 


CaSOLIMZ IKDEX 




0, 99 


0, 985 






Gasoline Yield (4) 




45 


44, 5 


43 


46 


Residue Yield (») 




12 


11,25 


13 


9,5 


ROK Clear Gasoline 




92,7 


92, 4 


92, 4 


92,4 



35 



45 
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TABLE 9.2 



"Sample" MCI obtained by densification 






MCI 


9B 


9D 




Fraction in Mixture 






0,889 


-0,276 
















A(cm-1) 


^(nm) 










4720 


2119 


0,001175932 


0,0012043 


0.00076051 


0,00081253 


4670 


2141 


0,017365599 


0.018971 


0,014737 


0,01 1 765 


4640 


2155 


0,017767782 


0,019076 


0.014281 


0,012239 


4615 


2167 


0,02411 8931 


0,025324 


0,019435 




4585 


2181 


0,019041012 


0,020612 


0,016021 


U,U 1 oii.0/ 


4485 


2230 


0.01 1987639 


0,011885 


0,010242 


0,010952 


4460 


2242 


0,017122848 


0,017086 


0,015362 


0,015912 


4385 


2281 


0,091829933 


0,089813 


0,09165 


0,096098 


4332 


2308 


0,1300549 


0,12812 


0,12815 


0.13281 


4305 


2323 


0,11911918 


0,11606 


0,1205 


0.12682 


4260 


2347 


0,09859607 


0,09716 


0,099018 


0,10194 


4210 


2375 


0,064520009 


0,063284 


0,066654 


0,068709 


4170 


2398 


0,061986802 


0,060439 


0,067026 


0,068963 


4135 


2418 


0,060186195 


0,058485 


0.065184 


0,06749 


4105 


2436 


0,062468569 


0,061975 


0,066333 


0,066187 


4060 


2463 


0.110237982 


0.11622 


0,11289 


0.098096 


4040 


2475 


0.093415832 


0,094281 


0,091751 


0.090004 














MONOO 


1 88,4 1 88.3 


87.2 1 87,6 
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TABLE 9.3 



"Sample" MCI obtained by densification 






MC2 


9B 


9C 


9D 


Fraction in Mixture 






1.162 


0,24 


-0.402 














X (cm-1) 


k (nm) 










4720 


2119 


0.001178718 


0.0012043 


0,00035436 


0.00076051 


4670 


2141 


0.017275892 


0.018971 


0,0048161 


0.014737 


4640 


2155 


0.01792139 


0.019076 


0,0062335 


0.014281 


4615 


2167 


0,023747506 


0.025324 


0,0088912 


0.019435 


4585 


2181 


0,019544078 


0,020612 


0,0084724 


0,016021 


4485 


2230 


0,011946806 


0,011885 


0.0093905 


0,010242 


4460 


2242 


0,017153608 


0,017086 


0,01448 


0,015362 


4385 


2281 


0,091428926 


0,089813 


0,099623 


0,09165 


4332 


2308 


0,13153274 


0,12812 


0,14239 


0,12815 


4305 


2323 


0,11829752 


0,11606 


0,13282 


0,1205 


4260 


2347 


0.098798684 


0,09716 


0.1071 


0,099018 


4210 


2375 


0,06434942 


0,063284 


0,073368 


0,066654 






0,061427506 


0,060439 


0,075591 


0,067026 


4135 


2418 


0,059908962 


0,058485 


0,075639 


0,065184 


4105 


2436 


0,062810524 


0,061975 


0,072756 


0,066333 


4060 


2463 


0,11139522 


0,11622 


0.090539 


0,11289 


4040 


2475 


0,09128022 


0,094281 


0.07754 


0,091751 














MONOO 




88.2 


88,3 


86 


87,2 
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TABLE 9.4 



Determination of MON of reformate based on "Samples" generated 






Weight 


Reformate 9A 


MCI 


MC2 








Measured 


Estimated 






Proximity Index 








0.000009271 


0.000010043 


0,000013457 
















A.(cm-1) 


k (nm) 












4720 


2119 




0,0010981 


0,0011773 


0,0011759 


0,0011787 


4670 


2141 




0,017744 


0,017321 


0,017366 


0,017276 


4640 


2155 




0,018144 


0,017845 


0,017768 


0,017922 


4615 


2167 




0,024297 


0,0239335 


0,024119 


0.023748 


4585 


2181 




0,020515 


0,0192925 


0,019041 


0.019544 


4485 


2230 




0,012619 


0,011968 


0,011988 


0.011948 


4460 


2242 




0,018197 


0,0171385 


0,017123 


0.017154 


4385 


2281 




0,092064 


0,0916295 


0,09183 


0.091429 


4332 


2308 




0,12886 


0,130795 


0,13006 


0,13153 


4305 


2323 




0,11882 


0.11871 


0,11912 


0.1183 


4260 


2347 




0,098322 


0,098697 


0,098596 


0.098798 


4210 


2375 




0,064577 


0,064435 


0,06452 


0,06435 


4170 


2398 




0,061405 


0,0617075 


0,061987 


0,061428 


4135 


2418 




0,059296 


0,0600475 


0,060186 


0,059909 


4105 


2436 




0,06198 


0,06264 


0,062469 


0.062811 


4060 


2463 




0,11037 


0,110815 


0,11024 


0,11139 


4040 


2475 




0,091698 


0,092348 


0,093416 


0,09128 
















MONOO 




0 


88 


88.3 


88.4 


88.2 
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TABLE 9.5 



15 



Recognition of new sample by auto adjustment 






Weight 


Measured 


Estimated 


9A 


Proximity Index 








0.000012235 


0.000012235 














A(cm-1) 


A(nm) 










4720 


2119 




0,0010702 


0,0010981 


0,0010981 


4670 


2141 




0,0171 


0,017744 


0,017744 


4640 


2155 




0.017768 


0,018144 


0,018144 


4615 


2167 




0,024103 


0,024297 


0,024297 


4585 


2181 




0.020269 


0,020515 


0,020515 


4485 


2230 




0.012224 


0,012619 


0,012619 


4460 


2242 




0.018338 


0,018197 


0,018197 


4385 


2281 




0,091998 


0,092064 


0,092064 


4332 


2308 




0,1306 


0,12886 


0,12886 


4305 


2323 




0,11841 


0,11882 


0,11882 


4260 


2347 




0,098802 


0,098322 


0.098322 


4210 


2375 




0,06262 


0,064577 


0,064577 


4170 


2398 




0,060234 


0,061405 


0,061405 


4135 


2418 




0,059762 


0,059296 


0,059296 


4105 


2436 




0,062527 


0,06198 


0,06198 


4060 


2463 




0,11151 


0,11037 


0.11037 


4040 


2475 




0,092677 


0,091698 


0,091698 














MONOO 






87.9 


88 


88 



40 



Claims 



1 . A method of determining or predicting a value of a property of a material X or a property of a product of a process 
45 from said material or yield of said process, which method comprises measuring the absorption Dj^ of said material 

at more than one wavelength in the region 600-2600 nm, comparing the said absorptions or a derivative thereof 
with absorptions Dim or derivatives thereof at the same wavelength for a number of standards S in a bank for which 
the said property or yield P is known, and choosing from the bank at least one standard with property P.^ said 
standard having the smallest average value of the absolute difference at each wavelength i between the absorption 
so D|X (or derivative thereof) for the material and the absorption Djm (or derivative thereof) for the standard S^, to obtain 
Fx. With averaging of said properties or yields P^ when more than one standard is chosen. 

2. A method according to claim 1 wherein the standard chosen for the property or yield wanted is such that in 
relation to the unknown material X and each chosen standard the following functions is met 

55 4 

< experimental error 
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wherein is the proximity index and is defined by i ^(xm) = i:(0 j^-D |^) ^ and the experimental error is in deter- 
mining said property or yield in the standard. 

3. A method according to claim 2 wherein the proximity index is less than the minimal index i^, which has been deter- 
mined from the standards Sa, Sb. So by (a) calculating for each pair of standards Sa/Sb, Sa/Sc the value of 

i2(a,b) etc, (b) relating the values of i2{a.b) etc to the corresponding differences EP (ab) in properties Pa, Pb etc (c) 
calculating for each value L for which i2(ab) is sL, the average of the corresponding differences EPab, (d) calculating 
Minimal index from the value of minimal index i2(ab} where average EPab is the same as reproducibility standard 
for the property. 

4. A method according to any one of the preceding claims, wherein the properties of synthetic standards, which are 
mixtures, and their spectra for consideration for possible choice for are estimated from existing standards in the 
bank for which, in respect of each existing standard for use in said mixture equation (4) and (5) are met, 

MinCj-uj s Cjg s MaxCj + Uj (4) 

and L Cj = 1 (5) 

wherein Cj is fraction of component j in the sample i 
Min Cj is the minimum of j in the samples for which the method is to be used 
Max Cj is the maximum of j in the samples for which the method is to be used 
uj is between l .0 and 0.05. 

5. A method according to any one of claims 2-4 wherein properties of Standards and spectra for consideration for 
possible choice are estimated by interpolation from measured properties of Standards and spectra for which the 
proximity index with respect to the unknown X is not more than 10 times the Minimal Index. 

6. A method according to claim 4 or 5 wherein at least one of (i) the estimated Standards and the corresponding 
spectra, and (ii) the property P^ of the unknown material and its spectrum, are added to the bank. 

7. A method according to any one of the preceding claims wherein the property is a physicochemical property of 

material X 

8. A method according to any one of claims 1 -6 wherein the property is a physicochemical property or yield of a product 
of a process to which at least one material X is a feed. 

9. A method for adding an extra synthetic standard to a bank of known standards, each of which relates at least one 
absorption in the 600-2600nm region (or derivative thereof) of a known material to a known property related to that 
material, which property is of said material which is a teed to a process, or is of a product from said process or yield 
of said process, which method comprises choosing from the bank at least 2 standards for which equations 4 and 5 
(as defined in claim 4) are met, considering mixing the chosen standards in at least one proportion to produce at 
least one mixture for use in a synthetic standard, and estimating the spectrum of said mixture according to equation 6 

S^^i = !:GjjXSj (6) 

where Sj is the spectrum in the mixture of component j in the calibration matrix, 
and estimating a property of said mixture according to equation 7 

PMi = i:CijXPj (7) 

where Pj is the property of component j, 

and then adding the spectrum and property of each "mixture" to the bank , and using them in at least one model 
involving a correlation/regression approach to relate NIR spectra to at least one property 

10. A method according to any one of claims 1-9 wherein said process is a hydrocarbon conversion or separation 
process, preferably a reforming or catalytic cracking or hydrotreatment, or distillation or blending. 

11. A . '->ethod according to claim 10 wherein said process is a distillation to give at least 1 distillation product and a 
residue and the properties/yields are obtained in respect of said product and/or residue. 
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1 2. A method according to any one of claims 1-11. wherein said property is in respect of a motor fuel and Is at least 
one of an Octane Number, vapour pressure, volatility, percentage distilled at 70 and at 1 0CC. gum content in 
mg/lOOml and content of sulphur, benzene, oxygen, methanol, ethanol, ter-butyl alcohol, methyl tert. butyl ether. 

TAME. 

5 

13. A method according to claim 12 wherein the property is in respect of a blend comprising gasoline the spectra are 
measured on feeds to said blending, and by calculation the blend index obtained as a linear or non linear function. 

14. A method according to any one of claims 1 -1 1 wherein said property is in respect of gas oil and is at least one of 
TO cetane index, cetane number, percentage of sulphur, density at 1 S^C, flash point, cloud point, f ilterability and viscosity 

at 40°C. 

15. A method according to any one of claims 1-10 wherein said property is in respect of a crude oil and is at least one 
of density, percentage of sulphur, viscosity at 100°C, content of paraffin, distillation curve and residual carbon per- 
is centage (Conradson Carton). 

16. A method according to any one of claims 1-10 wherein said property is in respect of a feed to a reforming process 
and is at least one of percentages of saturated linear, isoparaff ins, napthenics. and aromatics and density. 

20 M. k method according to any one of claims 1-10 wherein said property is in respect of a feed to a fluid catalytic cracking 
unit and is at least one of the density, the weight percentage of sulphur, the aniline point, viscosity at 1 00°C, refractive 
index at 20°C or 60°C, 50% distillation point, molecular weight, percentage of aromatic carbon and the KUOP, crack- 
ability or cokability of the feed or yield of gas, gasoline, gas oil or residue. 

25 1 8. A method according to any one of claims 1-10 wherein said percentage is in respect of the teed to a hydrogenation 
unit and is at least one of percentages of linear paraffins, isoparaffins, napthenes, linear olefins, cyclic olefins, dienes, 
benzene, toluene, xylene, alkylbenzene, density, or yield of light cut, heavy cut, or raffinate or benzene. 

19. Apparatus suitable for use in the method of any one of the preceding claims which comprises an NIR spectrometer 
30 receiving at least one signal from a feed or product line in said process and being coupled to a computer to effect 
continuous measurement of the spectra of the feed and/or product and provide feed back control of the process. 
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